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On  Thermal  Flffectc  in  the  Theory  of  Shells 


A.  E.  Green^  and  P.  M.  Naghdi 


Abstract . This  paper  is  concerned  with  thermomechanics  of  thin  shells  by  a 
direct  approach  based  on  the  theory  of  a Cosserat  surface  comprising  a two- 
dimensional  surface  and  a single  director  attached  to  every  point  of  the 
surface.  In  almost  all  previous  developments  of  the  thermo-mechanical  theory  of 
shells  by  direct  approach,  only  one  temperature  field  has  been  admitted.  This 
allows  for  the  characterization  of  temperature  changes  along  some  reference 
surface,  such  as  the  middle  surface,  of  the  (three-dimensional ) shell-like 
body,  but  not  for  temperature  changes  along  the  shell  thickness.  A main 
purpose  of  the  present  study  is  to  incorporate  the  latter  effect  into  the 
theory;  and,  in  the  context  of  the  theory  of  a Cosserat  surface,  this  is 
achieved  by  a recent  approach  to  thermomechanics  (Green  and  Naghdi  1977)  which 
provides  a natural  way  of  introducing  two  (or  more)  temperature  fields  at  each 
material  point  of  the  surface.  Apart  from  full  discussion  of  thermomechanics 
of  shells  and  thermodynamical  restrictions  arising  from  the  second  law  of 
thermodynamics  for  shells,  attention  is  given  to  a discussion  of  symmetries 
(including  material  symmetries)  and  thermal  effects  in  the  nonlinear  theory  of 
elastic  shells  with  detailed  discxission  of  the  linear  theory  of  elastic  plates. 
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1 . Introduction 


This  paper  is  concerned  with  thermomechanics  of  thin  shells  by  a direct 

approach  based  on  the  theory  of  Cosserat  (or  directed)  surfaces . A Cosserat 

surface  is  a body  c comprising  a two-dimensional  surface  (embedded  in  a 

Euclidean  3-space)  and  a single  director  (i.e.,  deformable  vector)  attached 

* 

to  every  point  of  the  surface.  A comprehensive  account  of  the  thermodynamical 
theory  of  a Cosserat  surface  --  hereafter  designated  C --  and  its  application 
to  shell  theory,  together  with  an  historical  survey  and  a large  number  of 
relevant  references  is  contained  in  an  article  by  Naghdi  (1972).  For  clarity's 
sake,  we  may  recall  that  the  material  surface  of  c can  be  identified  with  a 
particular  reference  surface  (often  taken  to  be  an  interior  surface)  in  the 
three-dimensional  shell-like  body,  e.g. , the  middle  surface  of  the  shell  in 
some  fixed  reference  configuration;  the  director  at  each  point  is  regarded  as 
representing  the  material  filament  across  the  reference  surface;  and  the 
component  of  the  director  along  the  normal  of  the  reference  surface  can  be 
taken  as  a measure  of  shell  thickness. 

Tliroughout  our  previous  developments  of  the  thermo-mechanical  theory  of 
shells  by  direct  approach  (Green  et  al.  1965,  Green  and  Naghdi  .1970,  Naghdi  1972), 
only  one  temperature  field  has  been  admitted  and  this  allows  for  the  characteriza- 
tion of  temperature  changes  along  the  reference  surface  of  the  shell-like  body. 

Tome  indication  of  how  temperature  changes  across  the  reference  surface  of  the 
shell-like  body  could  be  dealt  with  has  been  given  in  the  papers  of  Naghdi  (1964) 
and  of  Green  and  Naghdi  (1970,1971/  by  using  three-dimensional  approximations. 

One  author,  Zhilin  (1976),  has  considered  two  temperature  fields  in  a direct 
theory  and  we  refer  again  to  this  paper  below. 

» 

The  body  C is  taken  to  model  some  of  the  properties  of  a three-dimensional 
body  of  shell-like  character.  VJhen  the  director  is  absent  it  reflects  the 
properties  of  a material  surface  which  can  be  the  bounding  surface  between 
two  different  bodies  or  a surface  in  free  space. 
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Recently  mechanical  and  thermo-mechanical  theories  of  a membrane,  regarded 
as  a two-dimensional  surface  in  a Euclidean  3-space,  have  been  published  by 
Gurtin  and  Murdoch  (l97^»197b)  and  by  Murdoch  (1976a ,b),  where  the  surface  is 
mainly  regarded  as  a bounding  surface  of  a body  or  as  an  interface  between 
ambient  media  in  which  diffusion  through  the  interface  is  neglected.  Although 
these  authors  have  placed  special  emphasis  on  residual  stress  and  surface 
tension  in  the  membrane  surface,  apart,  from  notation,  all  their  basic  equations 
are  special  cases  of  those  given  previously.  Murdoch  (1976a),  however,  modified 
the  entropy  inequality  used  in  earlier  work  in  order  to  account  for  ambient 
temperatures  in  the  surrounding  media  which  may  be  different  from  each  other 
and  from  the  temperature  of  the  interface,  but  with  the  temperatxire  being  every- 
where continuous.  In  a second  paper,  Murdoch  (1976b)  seeks  to  show  that  his 
inequality  is  in  line  with  three-dimensional  considerations  of  Green  and 
Naghdi  (1970),  but  the  present  authors  believe  that  the  new  inequality  given 
by  Murdoch  is  unsatisfactory,  especially  when  only  one  temperature  field  is 
allowed  for  the  interface.  It  appears  that  there  may  be  confusion  between 
boundary  values  of  temperature  and  entropy  with  those  field  quantities  which 
should  enter  any  entropy  inequality.  We  believe  that  the  type  of  problem  con- 
sidered by  Murdoch  can  only  be  discussed  satisfactorily  on  the  basis  of  a two- 
dimensional  model  for  the  interface  if  more  than  one  temperature  field  is 
admitted.  In  the  context  of  a direct  approach  based  on  a Cosserat  surface, 
Zhilin  (1976)  formulates  a theory  in  which  two  temperature  fields  are  admitted; 
and  he  postulates  two  entropy  inequalities- , but  only  one  energy  equation  for 
the  surface.  Pince  this  provides  only  one  field  equation  for  the  two  tempera- 
tures, Zhilin  arbitrarily  rewrites  the  energy  field  equation  as  a set  of  two 
differential  equations.  The  physical  basis  for  these  two  equations,  which  have 
not  been  obtained  from  any  clearly  stated  balance  laws,  is  obscure.  It  is  dif- 
ficult to  see  any  relation  between  the  work  of  Zhilin  and  that  discussed  in 
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the  present  paper 


Although  widespread  use  of  the  Claus ius-Duhem  inequalities  has  been  made 
in  three,  two  and  one-dimensional  continuum  thermodynamics , these  inequalities 
nave  been  subject  to  the  criticism  that  in  some  circumstances  they  do  not 
reflect  adequately  ideas  associated  with  the  Second  Law  of  Thermodynamics. 

Green  and  Naghdi  (1977)  have  developed  a new  approach  to  three-dimensional 
continuum  thermomechanics  which  is  independent  of  any  particular  mathematical 
expression  of  the  second  law  and  which  imposes  some  restrictions  on  the 
constitutive  assumptions  leading  to  a reduction  of  a number  of  independent 
response  functions  (or  functionals)  in  the  set  of  constitutive  assumptions. 

In  the  present  paper  the  same  approach  is  used  for  the  Cosserat  surface  and 
this  provides  a natural  way  of  introducing  two  (or  more)  temperature  fields t 
When  the  director  is  absent,  the  theory  reduces  to  that  of  a material  surface 
which  may  be  a surface  in  free  space  or  a material  surface  between  two  dif- 
ferent media.  The  contrast  between  the  present  theory  and  that  of  Murdoch 
(1976a)  is  illustrated  by  an  example  in  §10.  On  the  other  hand,  for  an  elastic 
material  with  one  director  and  only  one  temperature  field,  we  recover  all  the 
previous  two-dimensional  results.  In  addition,  when  we  admit  two  temperature 
fields,  results  for  an  elastic  plate  agree  with  those  found  previously  by 
Green  and  Naghdi  (1970)  from  three-dimensional  considerations. 

Specifically,  the  contents  of  the  paper  are  as  follows.  Section  2 contains 
a concise  summary  of  the  various  basic  results  of  the  purely  mechanical  theory 
of  a Cosserat  surface  with  a single  director.  With  reference  to  thermal 
properties,  in  §3  we  admit  at  each  material  point  of  the  surface  of  C a number 
of  different  two-dimensional  temperatures  and  different  two-dimensional  entropies, 
as  well  as  related  thermal  fields;  and,  in  parallel  with  two-dimensional 


+For  the  purely  mechanical  theory,  it  is  already  known  hew  to  extend  the  theory 
with  more  than  one  director;  see,  e.g.,  Green  and  Naghdi  (1976). 
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conservation  laws  for  balances  of  mass  and  momenta,  we  postulate  balances  of 
entropy.  Next,  we  recall  the  balance  of  energy  for  the  Cosserat  surface;  and, 
following  the  recent  approach  of  Green  and  Naghdi  (1977),  after  elimination  of 
the  assigned  fields  --  i.e.,  assigned  force,  assigned  director  force  and 
external  rates  of  supply  of  entropy  — regard  the  resulting  equation  as  an 
identity  to  be  satisfied  for  all  thermo-mechanical  processes.  In  §4,  we 
briefly  discuss  thermoelastic  theory  of  a Cosserat  surface  on  the  basis  of  the 
new  procedure  in  thermomechanics  (see  §3)  and  also  compare  the  results  with 
earlier  developments  (see  Green  et  aJL.  196b;  Naghdi  1972)  involving  only  a 
single  temperature. 

A new  inequality  representing  the  second  law  of  thermodynamics  for  shells 
based  on  the  present  authors'  earlier  work  (Green  and  Naghdi  1977, 1978),  along 
with  restrictions  on  heat  flux  vectors  and  the  specific  internal  energy  are 
obtained  in  $5  and  §8,  respectively,  while  §§6-7  contain  a discussion  of 
relevant  results  for  shells  obtained  from  the  three-dimensional  theory.  The 
last  two  sections  (§§9-10)  are  devoted  to  a discussion  of  symmetries  (including 
material  symmetries)  for  shells  and  the  linear  thermoelastic  theory  of  isotropic 
plates.  The  developments  in  §§9-10  supplement  our  earlier  results  by  direct 
approach  (Green  and  Naghdi  1970,1971)  for  thermoelastic  shells  in  the  presence 
of  a single  temperature. 

The  general  theory  given  here  is  immediately  available  for  problems  in 
which  the  effect  of  surface  tension  and  interfacial  energies  in  a membrane 
surface  are  important,  but  we  do  not  consider  such  problems  in  detail.  Further 
discussion  is  necessary  if  there  is  diffusion  across  an  interface,  both  in  the 
context  of  the  membrane  theory  and  the  theory  of  a Cosserat  surface. 


4. 
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2.  Summary  of  mechanical  theory. 

We  summarize  in  this  section  the  main  kinematics  and  the  basic  equations 
of  the  purely  mechanical  theory  of  a Cosserat  surface  C and  refer  the  reader 
to  Nnghdi  (1972)  for  details  and  additional  references  on  the  subject.  Let 
the  particles  of  the  material  surface  of  (•  be  identified  with  a system  of 
convected  coordinates  0a  (or  =1,2)  and  let  the  surface  of  C in  the  present 
configuration  at  time  t,  hereafter  referred  to  as  J,  occupy  a two-dimensional 
region  of  space  R bounded  by  a closed  curve  3R.  Similarly,  in  the  present 
configuration,  an  arbitrary  material  surface  of  C occupies  a portion  of  the 
two-dimensional  region  R,  which  we  denote  by  P (cr)  bounded  by  a closed 
curve  3P.  Let  r and  d --  each  8 function  of  0°  and  t --  denote,  respectively, 
the  position  vector  of  a typical  point  of  J relative  to  a fixed  origin  and 
the  director  at  r.  Then,  the  base  vectors  along  the  0a-curves  on  J are 

/■w 

defined  by 


a = a (00,t)  = ar/de®  , 

''Of  ''■Of  ^ 


(2.1) 


and  we  denote  the  unit  normal  to  J by  a0(ea,t).  Also 


or  .a  or8  or  B 

a « = a • > 8 • a„  = 6,  , a ^ = a -a  , 

aB  ~a  ~8  ~ ~8  B ~ 


a = det  a . , [a.a0a_]  > 0 , a^a.  = a.  xa  , 

crp  /s<- 


where  6®  is  the  Kronecker  delta  and  aa  are  contravariant  base  vectors, 

8 

A motion  of  the  Cosserat  surface  is  defined  by 


(2.2) 


r = r(e®,t)  , d = d(0a,t) 


(2.3) 


5- 


r-  IP  I 


and  we  assume  that  d is  nowhere  tangent  to  Jy  so  that  a.,  ♦ d/0.  The  velocity 
and  the  director  velocity  vectors  are  given  by 


v = r , w = d , 


(2.4) 


where  a superposed  dot  denotes  differentiation  with  respect  to  t,  holding  d01 
fixed.  Throughout  this  paper  we  use  standard  vector  and  tensor  notations  (see 
Maghdi  1972).  ^Sreek  indices  take  the  values  1,2  and  the  usual  summation  con- 
vention over  a Greek  superscript  and  a subscript  is  employed. 

Consider  now  a reference  configuration,  which  we  take  to  be  the  initial 
configuration,  of  the  Cosserat  surface  C*  bet  the  reference  surface  in  this 

eonf igurat ion  be  referred  to  by  S with  K as  its  position  vector;  let  A and  A, 

~ ~o 

denote,  respectively,  the  base  vectors  along  the  0a-curves  on  S and  the  unit 
normal  to  J;  and  let  D be  the  reference  director  at  R.  Then, 


R « R(9*)  - r(Qa,0)  , D - D(6®)  - d(8“,0)  , A = dR/d9° 


(2.5) 


^ a ^ ^ a®  • a s s a®  • 

-or  * ~ 0*  ~ ~ ’ 

A * drt  A0*  • >0  • 


(2.6) 


We  assume  that  the  kinetic  energy  per  unit  area  of  the  surface  of  C in  the 
present  configuration  is  given  by 


A.  p 

T * ^ p(v.v+2k  v . w + k w • w ) , 


(2.7'> 


where  p*  Pi0Y,t)  is  the  mass  per  unit  area  of  J and  the  inertia  coefficients 
12  v 

k ,k  ‘ are  functions  of  9T  and  independent  of  t.  We  define  momenta  corresponding 


to  v and  w as 


. ■*>  * 


= e(v,+  k1*!  , = p(k*v  ♦ k*  w ) , (2.8) 

respectively,  per  unit  area  of  J. 

With  reference  to  the  present  configuration  at  time  t,  for  each  part  p of 
J , w postulate  the  equations  of  mass  conservation,  momentum,  director  momentum 
and  moment  of  momentum  as  follows: 


d r , l \ , 

— , piv  ♦ k w)dj 


pfio  * 


Nds 


P 


dP 


— - ' pvk^v  + k"w 'do  [ (pt  - m)do  f Mds 
iU  Jp  ~ ~ ~ ~ 


(2.0) 


(2.10) 


(2.11) 


_d_ 

dt 


r 

Jp 


p{r  x (v 


+ k1W) 


+ d x (k 


1 

V 


+ 


k'  w ' ]dj 


r - r 

I pir  x f + d x l 'do  + ' (r  x N ♦ d x M)ds 
JP J^P 


(2.1* 


o \ 


In  (2.10)  to  (2.12',  N - N(0\t  : \j)  is  the  force  vector  and  M M(d\t  ; u)  the 
director  force  vector^  at  the  curve  dp,  where  v is  the  outward  unit  normal  to 
dP  and 


or 


v a 


(2.13) 


The  vector  field  m in  (2.11' is  the  internal  director  force  per  unit  area  of  j 
and  it  makes  no  contribution  to  the  moment  of  momentum  equation.  Also 
£=£(0Y,t'  is  the  assigned  force  vector  and  L l(  0V,O  is  the  assigned 
director  force  vector  per  unit  mass.  These  include  contributions  which  model 
the  action  of  forces  over  the  major  surfaces  of  a shell. 


The  vector  N,  which  represents  force  per  unit  length,  has  the  dimension  [MT  ] 
where  [ M ] and  [T]  stand  for  the  physical  dimensions  of  mass  and  time.  If  the 
director  d has  the  dimension  of  length  then  k^,k‘  are  dimensionless  and  M has 
the  same  dimension  as  N.  If  d is  dimensionless  then  M has  the  dimension  t MLT”*- ] 
where  [ L]  is  the  physical  dimension  of  length.  In  the  latter  case,  M is  some- 
times called  a director  couple. 
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7. 


Under  suitable  continuity  assumptions  the  curve  force  vector  N and  the 


director  force  vector  M can  be  expressed  as 


N = N°v  , M ■ M°v  , 

Q[ 


(2.14) 


where  N0(,M0f  transform  ns  contravariant  surface  vectors, 
corresponding  to  (2.9)  to  (2.12)  are  then 

pa-'  * X(0Y)  or  p+  paa  • v =0 

~ ~,o 

(a^N°)  + Xf=X(v  + k1w)  , 

y Q(  <V  SW  <V 

(n'M®)  + \t  - mar  + X^k^v  + k^w ) 

^ yQl  *W  (V  »S# 

a xN°+dxm+d  xM°'=0  , 

'■'■O  ^ ^ ^ ~yQl  ~ 


The  local  field  equations 


(2.15) 

(2.16) 

(2.17) 

(2.18) 


where  a comma  denotes  partial  differentiation. 

In  the  absence  of  the  director,  the  field  equations  for  the  surface  are 
reduced  to  (2.15)  together  with 


+ Xf 

» Of  ~ 


Xv 


a x N®  = 0 


(2.19) 


These  are  one  form  of  membrane  equations  of  motion.  The  component  forms  of 
these  equations  were  noted  by  Green  et  al.  (1965,  §7)  and  were  derived  from  an 
energy  equation  together  with  invariance  conditions  under  superposed  rigid  body 
motions.  For  a more  general  derivation  using  invariance  conditions  and  a 
further  discussion  of  the  membrane  theory,  see  Naghdi  (1^72,  pp.  487-400  and 
pp.  546-547'.  If  an  interface  between  two  ambient  media  is  regarded  to  be  a 
membrane  surface,  then  the  membrane  equations  for  the  interface  must  be 
supplemented  of  course  by  appropriate  equations  for  each  media. 


? . '.hernial  properties.  Thermodynamical  theory  of  .-hells. 

In  most  existing  works  on  the  theory  of  a Cosserut  surface  only  one  temperature 
field  is  admitted  and  this  is  regarded  as  representing  the  temperature  variation 
ir.  some  reference  surface,  such  as  the  middle  surface,  of  the  sheil-like  body. 

-Iso,  the  effect  of  the  thermal  boundary  conditions  on  the  major  surfaces  of  the 
shell-like  body ^ are  incorporated  into  the  theory  through  the  external  surface  rate 
of  supply  of  heat.  The  variations  of  the  temperature  along  the  shell  tnickness 
have  not  been  model  led  so  far  by  a direct  approach  (within  the  scope  of  the  theory 
of  Cosserat  surfaces'',  although  some  indications  of  how  this  could  be  effected  is 
implicit  in  some  work  on  thermoelastic  shells  from  the  three-dimensional  equations 
by  the  present  authors  (Naghdi  1964,  Green  and  Naghdi  1970)  and  in  the  monograph 
by  Naghdi  (1972).  As  already  noted  in  §1,  because  of  the  new  approach  to  thermo- 
mechanics of  continua  introduced  recently  (Green  and  Naghdi  1977>1978),  it  is 
new  possible  to  account  in  a more  general  manner  for  the  thermal  properties  of  a 
shell-like  body  in  the  direct  formulation  of  the  theory  based  on  a Cosserat  surface. 

Thus , at  each  material  point  of  the  material  surface  of  C>  we  introduce  the 
scalar  fields*  0=0(0Y,t)  and  0jj=  0N(0Y,t),  (N  = 1,2,...  ,K) , representing  the 
effects  of  the  temperature  variation  in  a shell-like  body:  the  surface  tempera- 
ture 0,  which  we  require  to  be  positive,  represents  the  absolute  temperature  in 
the  reference  surface  of  the  shell-like  body,  while  the  scalars  0^,  account  for  the 

temperature  variations  along  the  thickness  of  the  shell.  In  addition  to  the 

* 


temperatures  0 and  0^,  we  admit  the  existence  of  external  rates  of  supply  of 

heat  r = r(0Y,t),  r =r  (0Y,t)  per  unit  mass  of  J and  external  rates  of  curve 
N IM 


supply  of  heat  -h,  -h^  per  unit  length  acting  across  the  boundary  dft*  Also,  we 

+The  terminology  of  major  surfaces  refers  to  the  upper  and  lower  surfaces  of 
the  body  separated  by  the  thickness  of  the  shell. 


* 


No  confusion  should  arise  from  the  use  of  the  symbol  0 in  the  uesignation  of 
the  temperature  fields  by  0, 0^ , ©2 , • ■ • , 0^  and  the  notation  0Y=(01,0i?''  for  the 
convected  coordinates. 


The  external  rates  of  supply  of  heat  r and  rN  include  contributions  corresponding  to 
heat  fluxes  on  the  major  surfaces  of  the  shell.  They  are  not  the  same  as  quantities 
defined  with  a similar  notation  in  Green  and  Naghdi  (1070'  or  in  Naghdi  (l^p^. 
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assume  the  existence  of  internal  curve  fluxes  of  heat  -h  = -h(g\t  ; v)» 

*■» 

-hj^  = -hN(  0Y,t  ; v)  across  each  curve  dP;  the  fields  h and  h^,,  called  heat  flxixes 
and  measured  per  unit  length  per  unit  timef  assume  the  values  h and  h^  on  dR, 
respectively.  The  total  external  rate  of  supply  of  heat  per  unit  mass  of  J 
is  defined  as 

K 

r + S r (3-D 

N=1  a 

Similarly,  the  total  external  rate  of  supply  of  heat  per  unit  length  per  unit 
time  across  dR  and  the  total  internal  curve  flux  of  heat  across  dp  per  unit 
length  per  unit  time  are  defined,  respectively,  by 

K K 

-h-  E h and  -h-  £ h . (3.2) 

N=1  N=1 

We  now  define  the  ratios  of  the  heat  supplies  r and  r^  to  temperatures  9 
and  8jj,  respectively,  as  s = s^0  ,t,  and  spj  = 9^>t ) and  call  these  the  external 

rates  of  supply  of  entropy  per  unit  mass  of  J.  Further  we  define  the  ratios  of 
h and  h to  the  temperatures  9 and  9N,  respectively,  as  the  external  rates  of 
curve  supplies  of  entropy  k and  k per  unit  length  of  dR;  and,  similarly,  we 
define  the  ratios  of  h and  h^  tc  the  temperatures  9 and  9^,  respectively,  as 
the  internal  curve  fluxes  of  entropy  k = k(9V,t  ; v)  and  k..=  k.,(9Y.t  ; v)  per  unit 

/-w  N IM  o- 

length  of  dP.  The  above  definitions  may  conveniently  be  summarized  by 


s=r/0  , 


SN  rN/0N  ’ k = h/9  » kN  ~ hl/®N  ’ 


(3.3) 


k’h/6  > W%  • 


We  require  that  the  fields  defined  by  (3.3)2  4 Q all  tend  to  finite 

limits  as  for  each  N = 1,2,...,K. 


The  sign  convention  for  the  internal  curve  fluxes  of  heat  are  such  that  these 
fields  -h  and  -hN  represent  fluxes  entering  P across  the  boundary  curve  dp. 


■ mm 
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In  addition  to  the  thermal  fields  already  introduced,  throughout  J we 
assume  the  existence  of  scalar  fields  TJ  = T}(  0^,  t ) and  Tlfj  = T1N(  9“,t ) , called 
specific  entropies  and  internal  rates  of  production  of  entropies 

and  »w=|N(0Y,t)  per  unit  mass  of  j.  The  contributions  of  these 
internal  rates  of  production  of  entropies  to  the  internal  rate  of  production  of 
heat  is 

K 

05 + £ eNsN  (3-4) 

N=1  1 

per  unit  mass. 

We  now  postulate  balances  of  entropy  for  every  material  surface  of  c occupying 

* 

a part  f3  in  the  present  configuration  and  write 


dt 


I p(s+§)d<7-  I kds  , H = | p T]  da  > 


L’<VVto-L>ds  = ■ “k-Jp’In'W  • 


(3-5) 

(3-6) 


By  usual  procedures,  it  can  be  shown  from  (3*5)  and  (3*6)  that  k and  are 
linear  functions  of  v,  i.e., 


k = p • v = p v 


kN  = hi  • a = PN  Va 


a 

p = p a 
— ~a 


PSSa 


(3-7) 


where  PjPjj  are  called  entropy  flux  vectors . Then,  from  (3*3)^  g and  (3-7)^  ^ 


h = ®E  * ~ » hN  = ' 't  (no  sum  on  N) 


(3-8) 


and  we  may  define  heat  flux  vectors  q, 


by 


% = ep  > Sn  = 0^  (no  sum  on  n) 


h = q • v 


hN  = 


Sn 


V 


(3-9) 


A motivation  for  postulating  equations  (3.5)  and  (3.6)  for  balances  of  entropy 
is  provided  by  consideration  of  derivations  from  three-dimensional  equations 
in  §7. 


11. 


Under  suitable  continuity  assumptions  and  with  the  use  of  (3.7),  the  field 
equations  resulting  from  (3*5)  and  (3.6)  are 


pT]  = p(s  + %)  - div  p , 

(3-10) 


where 


div  p = a'*(aV)  „ » div  = a‘^(a->“) 

~ ,ar  -J'J  N ,a 


We  now  introduce  the  first  law  of  thermodynamics  or  the  balance  of  energy 
for  the  Cosserat  surface  C.  This  states  that  heat  and  mechanical  energy  are 
equivalent  and  that  together  they  are  conserved  for  every  material  surface  of  C- 
Thus,  with  reference  to  the  present  configuration,  the  balance  of  energy  may  be 
stated  in  the  form^ 

j p 12 

— f 1 (v  . v + 2k  v • w + k w . w)  + e }pda 
at.  J p **  ^ ~ ~ ~ ~ 

r , K 

= (r+  E r + f . v + l • w)pdo 

**  p N=1  w **  **  **  ~ 

i*  ^ 

+ (N-v  + M.w-h-  E hjds  , (3.121 

J df>  ~ ~ ~ ~ N=1  N 

where  c = tl0Y,t)  is  the  internal  energy  per  unit  mass  of  J.  With  the  help  of 
(2.12)  to  (2.15),  (3.9)  and  (3.10)  and  under  suitable  continuity  assumptions, 
the  field  equation  resulting  from  (3.12)  is 

K K K 

- p(«-6V  z 0Ny-»ies<  r W-E- ft-  J = &,  ■ & 

N=1  N=1  N=1 

tN^'V  4m-w)M®'W  =0,  (3.13) 

***  ^ ~ ~ '**yQt 


where  the  temperature  gradients  g and  g are  defined  by 


^fee  Green  and  Naghdi  (197?)  for  further  remarks  on  this  in  the  context  of 
the  three-dimensional  theory. 


L 
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l) 


g grad  e 9 , g..  grad  9.,  9..  sa  . 

~ ,ar~  -M  N N,ai~ 

Introducing  the  He  Imho  1 tv.  free  energy  ^ ^r(9\t)  per  unit  mass  of  J by 


♦ C * 0T)  • E 9n1]n  , 

N 1 

the  energy  equation  it. 1*0  may  be  written  In  the  alternative  form 


(<.  l4) 


(3.1*>) 


-p^  ' ’19  ' E * Pv05  ' X &j^N  ~ 


+ Na  • v t m • w + M01  • w 0 

~ ~,Ur  ~ ~ ~ ''.or 


i t.H') 


For  a given  I'osserat  surface  having  a reference  density  ( 0a  ^ , the  fie  ld 
equations  obtained  from  the  integral  form  of  the  conservation  laws,  involve  a 
set  of  *' K * L *>  functions.  These  consist  of  the  deformat  ion  functions  r,d  and  the 

IN*  |N* 

temperatures  0.0N»  i.e. 


(?. 17) 


the  various  mechanical  and  the  thermal  fields,  namely 


(3.18) 


t£‘£’S,SN^ 


^ t.  I'll 


We  assvune  that  the  fields  (1.18)  are  specified  by  constitutive  equat  ions,  which 
may  depend  on  the  variables  (3.17),  their  space  and  time  derivatives,  as  well  a 
the  whole  history  of  deformation  and  temperature.  We  t lu-n  adopt  the  following 
procedure  in  utilising  the  conservation  laws: 


The  mass  density  p is  not  included  in  (3.18)  and  (3.1*0  since,  given  (i.iy), 
p can  be  calculated  from  (8.1*0. 
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(1)  The  field  equations  are  assumed  to  hold  for  arbitrary  choice  of  the 
functions  (3.17)  including,  if  required,  an  arbitrary  choice  of  the  space  and 


i 

M 

s \ 

|| 

* 

* 


time  derivatives  of  these  functions; 

(2)  The  fields  (3.18)  are  calculated  from  their  respective  constitutive 
equations ; 

(3)  The  values  of  the  variables  (3*19)  can  then  be  found  from  the  balances 
of  momenta  (2.1b)  and  (2.17)  and  balances  of  entropy  (3. 10); 

(4)  The  equation  (2.18)  resulting  from  the  balance  of  moment  of  momentum, 
and  the  equation  (3.16)  resulting  from  the  energy  equation,  will  be  regarded  as 
identities  for  every  choice  of  (3. 17).  This  will  place  restrictions  on  the 
constitutive  equations. 

We  note  that  the  quantities  5; , , T) ^ ^ may  be  arbitrary  to  the  extent  of 
• . K 

additive  functions  f , f ^ , f , f ^ , -0f  - £ O^f^,  respectively,  where  f,f^  are 

arbitrary  functions  of  the  variables  (3.17),  their  space  and  time  derivatives 
and  functionals  of  their  histories.  The  additive  functions  have  the  property 
that  they  make  no  contribution  to  the  differential  equations  for  r,d,0,0^  and 
the  boundary  and  initial  conditions.  They  also  make  no  contribution  to  the 
energy  identity  (3.1b)  and  no  contribution  to  the  internal  energy  s.  We  remove 
this  arbitrariness  by  setting* 

f fi0Y)  , fN  = fN(0Y)  , f = 0 , fN  « 0 . (3.20) 

Then,  the  functions  5,5^  are  determined  uniquely  and  T|,T^  are  only  arbitrary  to  the 

v A b 

extent  of  additive  functions  of  0,  independent  of  t.  The  functions  f,f^  in  (3.20)  can 
thus  be  determined  by  specifying  values  for  in  some  reference  state. 

Oo  far  no  mention  has  been  made  of  restrictions  on  constitutive  equations 
which  may  arise  from  some  form  of  second  law  of  thermodynamics , usually 

* For  a more  elaborate  parallel  discussion  in  the  context  of  the  three-dimensional 
theory,  see  Green  and  Naghdi  (ll}77,§2). 
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Interpret ed  In  terms  or  an  "entropy  Inequality,"  Before  considering  this  tout 

In  order  to  gain  pome  Insight  Into  the  nature  of  our  procedure  described  above, 

we  Study  tn,  t tie  next  section  the  relatively  simple  cnee  of  an  elastic  shell. 

For  later  use,  we  record  the  expressions  for  the  external  work  and  the 

external  heat  supplied  tv  any  part  s’  op  the  surface  ,/  during  the  time  Interval 

t,  * t * t . First,  however,  guided  by  the  results  of  $!*  we  observe  that  in  the 

case  of  an  elastic  Oosserat  surface  the  response  functions  .«  depend  only 

on  the  vectors  a ,d,d  and  the  temperatures  0,0.,  IN  l,?,...,iO  and  are 
^ ,y  N 

Independent  of  their  rates  and  the  temperature  gradients  g,g„.  i’uch  an  elastic 

m M\ 

material  will  be  regarded  as  nondl es i pat Ive  In  a sense  that  will  be  made  precise 
later;  and  In  conjunction  with  an  expression  for  the  external  mechanical  work 
supplied  to  any  part  P,  will  be  used  as  a basis  for  establishing  In  §'<  an 
Inequality  representing  the  second  law  of  thermodynamics  for  dissipative 
materials.  Keeping  this  background  in  mind,  wo  assume  that  the  constitutive 

response  functions  for  c,T]  include  also  dependence  on  the  list  of  variables 

* • • % % 

a « d , d -.9,0,  ,g,g_,  and  their  higher  space  and  time  derivatives  and  refer  to 

Ml  ^ A ** 

this  list  collectively  as  \/.  Further , let  c’.V  denote  the  respective  values  of 
c,T|when  the  list  V is  put  equal  to  cero  lit  the  response  functions.  Tints,  for 
example, 


C - C (a  ,d,d  ,0,0M,O  % 

mi ,y  N 


c'la  ,d,d  ,e,e'  # (•  »d ,d  ,0,0 N,o>  , 


O-n) 


V ( % d « d « 0 « * K % K*.  % • • • ^ 
^ y N ^ mNI 


where  the  dots  In  (U2l)  refer  to  the  higher  space  ami  time  derivatives  of 
• • • • • 

a^.d.d  p*  o..  .g  Then,  with  the  help  of  (2.1‘0  tv'  (2.1'(M  aitd  the  Integral 

of  (t.12)  with  respect  to  time,  we  obtain 


r 


Ei 


- 


ID  - eternal  mechanical  work  supplied  to  a part  p of  the  shell  during  the 
time  interval  t^  St  St,, 

* f2[f  (f  • v + 4 • w)p  do  + f (N  • v + M . w)ds]dt 
Jt  Jp  ~ ~ ~ ~ J dP  ~ **'  ~ ~ 


AK  > AE  ♦ Id  + Hi, 


(3.22) 


and 


M * External  heat  supplied  to  a part  p of  the  shell  during  the  time 
interval  tj  I t I t0 

pip  p ^ l»  K 

* [ 1 (r  + £ r )p  do  - I (q  + £ n ) . v dsjdt 

Jt1  V N-l  N 'dp  * N-l  '8N  ~ 


- ID- 


where 


(3-23) 

(3.24) 

(3-25) 


5"J"Jpl>(sf|'  + 1(fl9NVd<’dt  • *»  ■ • 

K and  E stand  for  the  kinetic  and  internal  energies  defined  by 
K » j k(v2  + 2k1v  • w + ktw2)p  da  , E = f p c do  , 

Jp  ~ ~ ~ ~ Jp 

respectively,  and  where  the  prefix  A denotes  the  difference  operation  on  functions 
and  fields  during  the  time  interval  [t-^tg],  e.g.,  AE  = E(t^)  - E(t^ ) . Also,  w in 
(3.24)2  is  given  by 

K 

pw=>  -p(c  - e ')  - p(*'  + Tl'0+  £ Tl'8.,)  + N®  • v +m.w  + Ma.v 

N«=l  N ~ ~>a  ~ ~ ~ ~*cr 

K K 

+ P(8?4  / °N5N)+E‘&  + h£  * (3.26) 

N=*l  N*=l 

*' " Ism  • ^.2?) 

N-l 

The  foregoing  discussion  of  thermodynamics  of  shells  includes,  of  course, 
the  thermodynamics  of  a membrane  surface  or  an  interface.  Results  of  this 
kind  can  be  recovered  by  deleting  the  terms  which  involve  the  director. 
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4.  Thermoelastic  theory  of  a Cosserat  surface. 

Previous  work  on  the  thermoelastic  theory  of  a Cosserat  surface  made  use 
of  a two  dimensional  Claus ius-Duhem  inequality  and  only  one  temperature  field 
was  considered,  which  corresponds  to  the  surface  temperature  9 of  the  present 
paper  (see  Naghdi  1972).  We  consider  now  constitutive  equations  for  a thermo- 
elastic  Cosserat  surface  which  admits  K+l  temperature  fields  and  we  examine 
the  restrictions  imposed  on  these  equations  by  the  procedure  described  at  the 
end  of  §3* 

We  assume  that  the  set  of  variables 


(4.1) 


are  functions  of  the  set 


, d , d , 9 , 0N , 8 , gN  , 

MY  ~ In  e-*  MM 


(4.2) 


as  well  as  the  reference  values 


A ,D,D  ,0  , 

~a  ~ ~,y 


(4.3) 


and  in  addition  may  depend  also  on  the  particle  9^.  In  the  set  of  reference 
values  (4.3),  0 is  the  constant  reference  value  of  9 and  we  have  assumed  the 
reference  values  of  9^  (N = 1,2, . . . ,K)  to  be  zero.  Postponing  the  restrictions 
to  be  imposed  by  the  invariance  requirements  under  superposed  rigid-body  motions 
and  recalling  the  procedure  outlined  in  §3>  the  energy  equation  (3.16)  is 
identically  satisfied  for  all  thermo-mechanical  processes  provided 


dt/3 g = 0 , = 0 


(4.4) 
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,®  ; P)  , 

(4.5) 

£ = p ’ if  = p ’ 

~o  ~,cr 

<^1 

CL 

u 

Bl 

(4.6) 

A A 

-n  _ _ M T,  = _ 

" ' de  ’ \ d9N 

> 

(4.7) 

K K 

p(e?\;,  Vi’t-ri,  e» 

II 

O 

(4.8) 

N=1  N=1 


where  as  indicated  in  (4.5)  the  function  y is  independent  of  the  temperature 

gradients  g and  g^.  Formulae  (4.6)  and  (4.?)^,  with  9^,  absent  from  (4.5), 

were  obtained  by  Naghdi  (1972,  Sec.  13)  with  the  help  of  the  Claus ius-Duhem 

inequality  and  he  also  showed  how  to  obtain  alternative  forms  with  the  help  of 

invariance  conditions  under  superposed  rigid  body  motions.  In  this  connection, 

it  may  be  recalled  that  in  the  discussion  of  some  forms  of  the  constitutive 

equations  for  an  elastic  Cosserat  surface,  Naghdi  (1972,  Sec.  11 ) at  first 

indicated  the  dependence  of  the  response  functions  on  the  properties  of  a 

physically  preferred  reference  state  and  material  inhomogeneity  through  the 

argument  9^  and  subsequently  specified  a more  explicit  dependence  on  the 

reference  values  such  as  (4.3)  in  order  to  obtain  other  forms  of  constitutive 

* 

equations  in  terms  of  relative  kinematic  measures.  Results  for  an  elastic 
membrane  or  an  interface  Regarded  as  a membrane)  follow  from  (4. 5)  to  (4.8) 
by  suppressing  the  directors  D,d  and  omitting  the  response  functions  for  11°, m. 

It  is  now  convenient  to  introduce  the  component  form  of  the  kinematic 
variables  d,d  relative  to  the  base  vectors  a.  or  a*,  and  the  component  form 
of  P,D  relative  to  A.  or  A . Thus,  we  write 


The  influence  of  the  reference  geometry  on  the  response  of  elastic  shells  has 
been  examined  also  by  Carroll  and  Naghdi  ^1972^  who  assumed  the  existence  of  a 
local  preferred  state  of  the  body  and  then  stipulated  that  the  influence  of  the 
reference  geometry,  as  in  (4.5),  occurs  through  the  values  of  the  constitutive 
variables  in  the  preferred  state. 
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r 


UHN 


__________ 


where 


n 

d a* 
1<^ 

= d1^  , 

d = X.  a1  = 
~,a  ia~ 

x'c£i 

9 

s 

(4.9) 

D = 

D A1 
i~ 

- DlA.  , 

D * A.  A1  = 
~,a  ia~ 

A1  A. 

• 0^1 

9 

• 

d . = d . a 

1 ~ M. 

> 

d1  - 

o1Jd  , 

J 

X.  = d • n, 

~,a 

, x1 

•a 

-sl%.  ■ 

♦ 

Di = £•  *d 

9 

D1  = 

.ijn 

A Dj  , 

A,  = D -A. 

~,  Q 

-aUv  - 

(4.10) 

ai3  = 8 

3i  = 

0 , 

a33  = ]_ 

, Ai3=A3i= 

0 , A 

33  . 1 

and 


= d i - b d 

, X.„  - d„  + b^d 

<*P  a|p  aP  3 

’ 3®  3, a a p ’ 

®P  ‘'arllP  ‘ ap1  3 

, A-  = D_  + B^D  . 

oa  3. a a P 

(4.11) 


In  the  above  formulae,  a single  vertical  bar  stands  for  covariant  differentiation 


using  Christoffel  symbols  formed  from  a a double  vertical  line  denotes 

up 


covariant  differentiation  using  Christoffel  symbols  formed  from  A „ and 

b or0 


b _ = a • n _ , B . = A_  • A 
ap  ~3  ~ar.  P aP  ~3  ~or.P 


(4.12) 


are  the  coefficients  of  the  second  fundamental  form  of  the  surfaces  J and  S, 
respectively. 

With  the  help  of  invariance  conditions  under  superposed  rigid  body  motions 


the  constitutive  equations  (4.5^  to  (4.7)  may  be  expressed  in  an  alternative 


form  which  wiLl  be  utilized  later  in  the  paper.  Under  such  motions  the  vectors 


a ,d,d  become  Q a , Q d,  Q d , where  Q is  a proper  orthogonal  function  of 

~or  ~ ~ ~o  ~ ~ ~ v ~ 

+ A 

the  time,  and  the  value  ^ of  the  response  function  ^ in  (4.5^  is  given  by 


'See  Nnghdi  (1972,  Sect.  13 ). 
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i ~ ♦(ft  ® > Q.  d,  Q d , 9,  8,.,  A , D,  D , 8 ; o**) 

~ ~o<  ~ ~ ~ ~»V  n ~o  ~ ~>V 

We  assume  that  the  value  of  the  free  energy  is  unaltered  by  superposed  rigid 
body  motions  so  that  - $ or 

(Q  b»  Q d,  Q d ,.••)  = \ji(a  ,d,d  »•••) 

~ ~or  ~ ~ ~ ~,y  ~ ~.v 


(4.13) 


~Q 


(4.14) 


for  all  proper  orthogonal  tensors  Q.  At  this  point  it  should  be  mentioned 
that  those  who  prefer  the  concept  of  objectivity  will  make  the  statement  (4.14) 
for  the  full  orthogonal  group. 

Assuming  first  that  (4.14)  holds  only  for  the  proper  orthogonal  group  for 

A 

Q then  from  Cauchy's  representation  theorem  \[i  may  be  expressed  as  a different 

function  of  the  inner  products  and  scalar  triple  products  of  a ,d,d  , namely 

~a ,y 


1 


a • a.  = 

-or  ~0 

a , a 

0(0  ~o 

d = d 
~ a 

’ 4 ' 

d = XQ  , 
~,a  0o( 

d = dXd, 

^ 1 

, d • d 

~ ~jQ ( 

d\ 

ia 

* ~,a‘ 

d a = X1  X.Q 
0 • a 10 

(4.15) 


and 

[2a2ei]  = a?d3  » [Sd2e£,a]  = a%a  ’ 

[a  dd  ] = d3e  XY_  + \3  e dy  . 

~a  ~ ~,  0 ya  *0  .0c yy 

A A 

In  view  of  (2.2)  , [a, a_a, ] = a2  > 0,  so  that  a2  can  be  expressed  as  a single- 
5 

valued  function  of  a a.  It  then  follows  from  (4.15)  and  (4.16)  that  may  be 

ap 

expressed  as  a single-valued  function  of  a^jd^X^,  the  remaining  contributions 

in  (4.16)  being  redundant.  Similarly,  if  we  make  the  other  choice  of  the  normal 

a^  in  (2.2)  so  that  [a.a.a,]  = - a2  , we  may  again  reduce  * to  be  a single-valued 
~3  ~o 

function  of  a „,d.,X.  . Moreover,  when  [A-,A„A,]>0  or<0,  we  may  express  A_  as 
a0  i iQr 

A 

a single-valued  function  of  A , and  replace  the  dependence  of  4 on  A ,D,D  by 

'■'■Of  '*JQt  P 


(4.16) 
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A ot’  ^ence  (4-5)  is  replaced  by  the  different  single-valued  function^ 

* = ^naB,di,Xicr,0,0N  5 £a,Di  ,Aka’®  ; ' 


(4.17) 


Let  N0(i,M“i  ,ni  denote,  respectively,  the  components  of  N®,M®  and  m referred  to 


the  base  vectors  a , i.e., 
'—1 


N®  = N®V  , M®  = M®V  , m = mV  . 

~ ml  *+*  ml  ~i 


(4.18) 


i 


Then,  with  the  help  of  (2.2)  and  (4.10),  it  follows  from  (4.5),  (4.6),  (4.7)  and 
(4.17)  that 


■l-’&  - ■ 

N'®0  = N®0  - m®d0  - M^®\0  = 0(r^~  + -^-)  , 

^aB  ** 

N®3  - m®d3  - M^®\_  + m3d®  + MY~  \a  = 0 , 

3y  *y 


(4.19) 


n n =.  Ji. 

^ de  ’ ‘n  aeN  * 

Apart  from  the  last  formula  in  (4.19),  the  above  results  are  equivalent  to  those 
given  by  Green,  Naghdi  and  Wainwright  (,1965).  The  moment  of  momentum  equation 
(2.18)  is  satisfied  by  (4.19),  1 • 

In  order  to  complete  the  discussion  of  invariance  under  superposed  rigid 
body  motions  it  remains  to  consider  p,p„,£,5„.  For  example,  recalling  (3.7).,, 

'>■'  MM  lM  jf 

we  have  the  constitutive  relation 

P = P ^ 

^ ~ MM  /sw 


subject  to  the  condition 


'The  response  function  ^ in  (4.17)  exhibits  explicitly  dependence  on  the  basis 
in  the  reference  configuration  of  the  surface  of  C-  Since  j is  a scalar- 

valued  function,  its  dependence  on  A will  be  through  A =A  • A . 

~or  oB  ~a  ~B 
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p <a  v a • sa,y'  »>  v a &>  %sn<  •••> 

■ 5‘i!i#«,Y’,’Va;  -1 

for  all  proper  orthogonal  tensors  Q.  By  a discussion  similar  to  that  used  for 
>jf , it  follows  that  pa  can  be  expressed  in  the  different  functional  form 

p“  = f”<‘Wdl>’W6>V  7Z’  -Z'-  WV01  ^ • 

00  00 

Similar  results  follow  for 

Although  we  adopt  the  representation  (4.17)  in  the  rest  of  this  paper,  we 

note  that  if  (4.13)  is  to  hold  for  the  full  orthogonal  group,  then  \|i  becomes  a 

function  of  the  inner  products  (4.15).  Moreover,  since  a • d = d /0,  we  may 

A ~ ~ 

reject  redundant  elements  in  (4.15),  and  reduce  ji  to  a function  of 

aa0,da,Xae’(d3)2’d3>k3a  ; VDi,Aka’0  5 

as  found  by  Naghdi  (1972,  Eq.  (13.36)).  On  the  other  hand,  the  component  d, 
(since  it  is  nonzero)  can  either  be  chosen  to  be  always  >0  or  always  <0.  In 
each  case  d^  may  be  expressed  as  a single-valued  function  of  (d^)  so  that  ^ as 
a function  of  the  variables  (4.20a)  can  be  expressed  in  the  form  (4.17)  and 
conversely.  This  means  that  there  is  no  difference  in  results  if  we  use  the 
orthogonal  group  instead  of  the  proper  orthogonal  group  for  our  invariance 
conditions.  Moreover,  there  are  different,  but  equivalent,  representations 
for 

Returning  to  our  main  objective,  with  the  help  of  (4.8)  and  (3.22)- 
(3.24),  for  an  elastic  Cosserat  surface  we  record  below  the  expressions 
representing  (i)  work  by  assigned  force  and  assigned  director  force  and  by  con- 
tact force  and  contact  director  force  acting  or.  any  part  p and  (ii)  supply  of 
energy  arising  from  the  total  external  rate  of  supply  of  heat  and  the  total 
curve  flux  of  heat  to  P,  both  over  a finite  time  interval  t St*t0: 

•4  cL 


i 


IP  = f '"[  f (f  ■ v*  ( • w)p  + [ (N  • v + M • w)ds  ]dt 

V J3p 


= AK  + AE ' + lc 


( 


and 


: f "L  [ (r  + E rN)p  da  - [ (q  + Z q^)  • v ds 
Jt1  Jp  N=1  H Jdp  ~ N=1  ~ 


]dt 


= lc 


where  to  and  K are  given  by  (3-24).,  and  (3-25),,  respectively,  and  E'  is  defined 


5.  The  second  law  of  thermodynamics . 

Previously  (Green  and  Naghdi  1977,1978),  in  the  context  of  the  three- 
dimensional  theory  we  have  discussed  the  nature  of  thermodynamic  irreversibility 
“rising  from  a mathematical  interpretation  of  a statement  of  the  second  law  of 
thermodynamics  that  "it  is  impossible  completely  to  reverse  a process  in  which 
energy  is  transformed  into  heat  by  friction."  Here  we  follow  the  same  procedure 
and  reconsider  a mathematical  interpretation  of  a second  law  appropriate  for 
a Cosserat  surface  which  admits  more  than  one  temperature  field.  Earlier  work 
or.  the  subject  made  use  of  a Claus ius-Puhem  inequality  when  only  one  temperature 
field  is  admitted  (Green,  Naghdi  and  Wainwright  l^o^,  Naghdi  1972  . 

For  the  sake  of  clarity,  we  first  dispose  of  some  definitions.  A state  of 
the  Cosserat  surface  C at  time  t,  regarded  as  representing  a thin  shell-like  body, 
is  described  by  the  position  vector  r and  the  director  d,  the  velocities  (2.4)  , 

~ '■*>'  1 yd 

the  temperatures  0 and  0^  (N = 1,2, . . . ,K)  throughout  the  surface  J of  c,  together 
with  the  constitutive  response  functions  for  the  fields (3. IS) • Once  the  response 
functions  are  given,  we  then  know  the  values  of  H and  in  (3*5)2  and  (3-6),,, 
as  well  as  K and  E in  (3*25).  A thermo-mechanical  process  or  simply  a process 

1 yd 

is  a time  sequence  of  states;  it  is  a continuous  oriented  curve  in  the  space  of 
states,  i.e.,  the  ( 0,  ®yj,aa|3’di  ,4^-spnce.  Thus,  a process  may  be  defined  by  a 
sequence  of  values  of 


(e,eN,r,d) 


(5-1) 


throughout  J in  the  time  interval  05t5c.  Similarly,  the  reverse  process  is  a 
process  defined  by  a sequence  of  values  of  (5*1)  throughout  J in  the  time 
interval  a^tS2a  subject  to  the  conditions 


e(t)  = 0(2a-t)  , eN(t)  = 0N(?a-t)  , 


r(t)  = r(2a-t)  , d(t)  = d,.(2o-t) 


(5*2) 
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Returning  to  our  main  objective,  we  observe  that  in  any  process  the  work 
done  by  the  external  mechanical  forces  acting  on  p is  positive  or  negative 
depending  on  whether  the  external  work  is  supplied  to,  or  is  withdrawn  from,  p. 

In  general,  some  of  the  work  done  results  in  a change  of  the  kinetic  and  internal 
energies  represented  by  the  first  two  terms  on  the  right-hand  side  of  (3.22)^, 
each  of  which  may  be  positive,  negative  or  zero.  Also,  part  of  the  work  done 
may  be  positive  with  a corresponding  extraction  from  P as  heat  or  negative  with 
a corresponding  absorption  of  heat  by  p.  We  note  that  in  the  case  of  an  elastic 
material,  the  different  contributions  to  lb  will  vary  in  sign  depending  on  the 
process  and  will  not  be  restricted  to  be  either  positive  or  negative  for  all 
processes.  Consider  any  smooth  process  in  the  time  interval  OStSo  and  its 
reverse  process  in  the  time  interval  a^tS2a.  If  the  process  is  reversed  in 
the  (9, 0^,a^,d^ ,\^)-space  in  such  a way  that  at  the  end  of  the  process  and 
its  reverse  process  the  shell  has  returned  to  its  original  state  with  A9  = 0, 
A9n=0,  Aa  =0,  = °>  ^ia=0j  Av  = o,  Aw  = 0 and,  hence,  Ae  = 0,  AT)=0,  AT]N  = 0, 

A?  = 0,  A|M  = 0,  ANa=o,  AM°=o,  Am  = o,  AP  = °,  AP„=°  an(i  AK  = 0,  AE=0,  then  all 
the  work  done  in  the  process  is  recovered  as  work  in  the  reverse  processt  This 
recovery  of  work  would  not  be  possible  if  in  every  arbitrary  process  part  of  it 
always  has  a positive  sign.  With  this  motivation  in  mind,  we  assume  that  for 
any  arbitrary  process  in  a dissipative  shell  only  part  of  the  work  done  is 
recoverable  as  work  in  a reverse  process,  the  rest  being  transformed  into  heat. 

We  therefore  assume  that  in  every  process  part  of  the  work  done  is  always 
nonnegative.  Then,  if  at  the  end  of  any  process  and  its  reverse  process  the 
shell  has  returned  to  the  same  state,  some  of  the  work  done  is  always  trans- 
formed into  heat.  Recalling  that  Ut,  = 0 in  (3*24)^  in  the  case  of  an  elastic 
shell  for  all  processes,  we  interpret  the  above  assumption  for  a dissipative 


'if  work  is  extracted  in  the  process,  then  it  is  absorbed  by  the  medium  in  the 
reverse  process. 
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\Q?  * 0 (5.3) 

for  all  parts  p and  all  processes,  where  Uv,  is  given  by  (3-24)o.  Since  ,t 
are  arbitrary  and  pw  has  already  been  assumed  to  he  continuous,  it  follows 
that 

K 

pw  = - pt«-«  ')  - p(i'  + il'e  f X!  T) ' 9 ^ + N • v + m • w n Ma  • w SO  (s.4) 

in  — '■',0'  ~ ~ ~ a 

for  all  thermo -mechanical  processes.  Also,  from  (3.23)  and  (5.4),  we  have 

p*op  . N 

»s'  p(eV  + s eMiU)da  dt  (5.5) 

Jt1Jp  u"" 

so  that  the  external  heat  supplied  to  a part  p of  the  surface  of  C if5  hounded 
above  in  any  process. 
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urnmary  of  results  from  three-dimensional  mechanical  theor 


Consider  :•  three-dimensional  body,  embedded  in  a Euclidean  3*sPuce,  and 


let  its  particles  be  identified  by  a converted  coordinate  0 


In-'iis  Iona  i t ody  in  the  present  configuration  at  time  t.  Then 


where  g.  and  g are  eovariant  and  contravariant  base  vectors,  respectively 


and  g are  covariant  and  contravariant  metric  tensors,  respectively,  and 


is  the  Kronecker  delta.  Also 


superposed  dot  denotes  material  time  deriva 


The  stress  vector  t across  a surface  in  the  present  configuration  whose 


\init  outward  normal  is  n is  given  by 


where  t are  the  contravariant  components  of  the  symmetric  stress  tensor.  We 


do  not  recall  here  the  consequences  of  the  conservation  laws  of  the  three 


dimensional  theory  since  they  will  not  be  used  explicitly  in  the  present  paper 


The  parametric  equation  9J = 0 defines  a surface  in  space  at  time  t,  which 


we  assume  to  be  smooth  and  non-intersecting.  Any  point  of  this  surface  is 


specified  by  the  position  vector  r,  relative  to  the  same  fixed  origin  to  which 


r is  referred,  where 


Let  the  boundary  of  the  three-dimensional  continuum  be 


surfaces 


For  convenience,  we  adopt  the  notation  for  r in  ( 2* . 3 ) ^ also  for  the 
This  permits  an  easy  identification  of  the  two  surfaces. 


T 


I- 


e3  = ^(e1,©2)  , e?  = p(e\e2)  , (O<o<p)  , 


(6.4) 


with  the  surface  8*0  lying  entirely  between  (6.4)  and  a material  surface 

l ><- 


f(eL,e2)  = o 


(6.5) 


which  is  such  that  Q'' = const,  are  closed  smooth  curves  on  the  surface  (6.5). 
Since  later  we  identify  the  surface  described  by  (6.3)  with  the  surface  (2 . 3 ) } , 


it  is  convenient  to  adopt  the  notations  for  a ,a  , etc.  defined  by  (2.1)  and 

~a  arp 


(2.2)  also  for  the  surface  (6.3). 


We  may  now  consider  a general  representation  for  r in  (6.1)^  as 

3 


polynomial  in  0 but  in  what  follows  we  restriction  attention  to  the 
approximation 


r = r + 03d  , d = d(0Qf,t) 


(6.6) 


in  the  bounded  region  ai0JS0.  Recalling  (6.1),  the  velocity  vector  is 
then  given  by 


* 3 

v = v + 0'  w , w = d 


(6.7) 


Given  the  approximation  (6.6^,  it  is  known  (see,  e.g.,  Naghdi  1972;  Secs.  11, 
12)  that  the  field  equations  of  the  forms  (2.15)  to  (2.18)  can  be  derived  from 
the  three-dimensional  field  equations  provided  we  identify  d in  (6.6)  with  (2.3), 
and  adopt  the  definitions 


i (**  * -^  3 N 

Pa'  = k ® j P , \k  * I p g‘  (©’)  d 0 , 


(0.8) 


N°tr'  = f Tad03  , Maa'-  = f T°03d03  , nwr  = f T3d03  , (6.9) 


h rp 

Of 


where  0 is  the  three-dimensional  mass  density.  Also  the  assigned  force  f and 


H 
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the  assigned  director  force  l are  related  to  the  three-dimensional  body  force 
* 

f per  unit  mass  and  to  the  effects  of  the  stress  vector  (6.2  )^  over  the 
boundary  surfaces  (6.4)  0 by 

X y£- 

= fP  pVg^de3  + [tg^i  U3)]  , +[tg2fa(e3)]  , 


where 


•j  ~ 

cr 


Xj  of ' ^ * 3 ~ 3 ’ 

d -a  p 0 -p 


u - f PVgVW  + [tgVt  (e3)]  . + [tgVf  (e3)]  , 

~ ~ “ e3=«  ~ 0 e3^ 


f (e3)  - [gU(^)2  + g22(^)2  + g33 

ae  ae 

+ 2(g12  ^ *l  _ gi3  . /3  ^}]h 

ae  ae"  ae  ae 

and  f (93)  is  obtained  from  (6.9)  by  replacing  or  by  {}. 


7.  Thermodynamical  results  from  throe-dimensional  theory. 

In  this  section  we  obtain  some  thermodynamical  results  for  a shell-like 
body  on  the  basis  of  the  recent  thermodynamical  theory  of  Green  and  Naghdi 
(1977).  Thus,  along  with  the  (three-dimensional)  temperature  field  * 

9 =9(0  >t)  > 0,  we  admit  the  existence  of  an  external  rate  of  supply 

— * » 

of  heat  -h  per  unit  area  acting  across  the  boundary  dR  of  a region  of  space 

» 

R occupied  by  the  body  in  the  present  configuration  at  time  t.  Also,  we  assume 

* # i 

the  existence  of  an  internal  surface  flux  of  heat  -h  =-h(0,t;n)  per  unit 

» * 
area  across  each  surface  dp  which  is  the  boundary  of  an  arbitrary  part  p of 

* # * * * i * 

ft  . We  define  the  ratio  of  the  heat  supply  r to  temperature  9 as  s = s (0  ,t) 

and  call  this  the  external  rate  of  supply  of  entropy  per  unit  mass.  Similarly 

— # » 

we  define  the  ratios  of  h and  h to  temperature,  respectively,  as  the  external 

— * * 

rate  of  surface  supply  of  entropy  k per  unit  area  of  aft  and  the  internal  surface 

* *,  i * 

flux  of  entropy  k =k  (9  ,t  ; n)  per  unit  area  of  dp  . Thus, 

* * * — * » * * # 
r = 0 s , h = 0 k , h = 0 k . (7.1) 

* * * i 

In  addition,  throughout  ft  we  assume  the  existence  of  a scalar  field  T)  - T)  (0  ,t) 

per  unit  mass,  called  the  specific  entropy  and  an  internal  rate  of  production  of 
* *,  i \ 

entropy  5 = 5 (9  ,t)  per  unit  mass.  The  contribution  of  the  latter  to  the  internal 

# * 

rate  of  production  of  heat  is  simply  0 Z,  per  unit  mass. 

We  recall  the  balance  of  entropy  in  the  form 

dr**!,*/**\r*  . . 

-rr  p T1  dv  = p (s  + I )dv  - 1 k da  (7.2) 

dt  v V "dp 

* 

for  every  material  volume  occupying  a part  p in  the  present  configuration.  It 

» 

follows  that  k is  linear  in  n,  i.e., 

**•**} 

k = p«n  , p = p g.  , (7.3) 

rw  rss  f\u. 

* ...,*** 
where  p is  the  entropy  flux  vector.  Then,  flora  (7.1)  and  (7*3) > h = 9 p • n 

*No  confusion  should  arise  from  the  notation  (0l,9‘",9^ ) for  the  convected 

coordinates  and  the  use  of  the  symbol  9 in  designation  of  the  (three-dimensional) 

temperature  field  9*1  and  the  surface  temperatures  0,91  >0O, . . . ,0  in  (7.17). 
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and  we  may  define  the  heat  flux  vector  q by 


* * * 
q = 0 p 


(7-4) 


With  the  h.'lp  of  (7.3)  5 the  field  equation  corresponding  to  ( 7.2 ) is 


* • * * , * * 


P Tj  = p (s  + £ ) - div  p , 


(7-5) 


where 


div  p = g~‘:d(p',kg?)/dek 


(7-o) 


Now  multiply  (7*5)  by  (03)N  and  integrate  over  an  arbitrary  part  P in  the 
present  configuration.  After  using  (7-3)  and  some  straightforward  manipulation, 
we  obtain 


df  * *,  3sN  p *.  * *s , 3xN  P / 3,N-1  *3 
jT  P T)  O'3)  dv  = I p (s  )(03)Wdv+  N(0J)H  p 3dv 

at.  .ip-  p*  jp 

- f ,.(e3)Vda  (N  s 1,2, ... ) . (7*7) 

'!dP 

Let  an  arbitrary  material  surface  0=0  occupy  a region  P at  time  t and 

* * 

let  dP  denote  the  closed  boundary  of  p.  Further,  let  dP^  refer  to  a part  of  dP 

specified  by  the  surface  (6.5)  so  that  dP  = dp  = dp  on  0-0,  and  let 

n 

dP*c  = dP  - dPn  stand  for  the  complement  of  dP^  in  dP  . Then,  for  a shell-like 

region  bounded  by  the  surfaces  (6.4)  and  (6.5) > from  (7*2)  and  (7-7)  with 

1 >2 

N=1,2,...,K,  we  can  derive  the  balance  equations  (3.5)  and  (3-6)  withoiit 
introducing  any  approximations^  provided  we  make  the  following  identifications: 

P T)  g^0J  » \\  = ! 

a * a 


P*  * * i 3 i'P  » * 1 3J  -s 

\T\  = J P T)  g?d0-  , \T1n  = j p T\  g5(03)  d0“  , 


(7.3) 


(•0  * # 1 Q 111  o » 1 o 

\s  = I p s g5d0  -[k  g?t(03)]  -[k  g5ffl(0-)]  , 

a a 03=o  P 03-p 


(7-9) 


^The  details  parallel  similar  developments  in  Naghdi  (1972,  Sec.  11 ). 
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8' =a 

(•P  * * 1 3 

x?  = j p 5 e d9'  > 

a 


X...  f V.*gW)W-(k*e*(,>3)Bt  (e3)]  , -[k*gi(e3)"t,(»3)]  , , (MO) 

(7-11) 
(7.12) 

p“.*  - f p""gW  , P“ai.  J0pV(e3)V  , (7.i3) 


»P  * * i. 


>S»*i*rAV(^)%3  , - n f*  pV(03)n-V  , 


(7.14) 


(7.15) 


k = p\  ’ kN  = pNva  ’ 
h = 9k  , hN  = 9^  . 

We  also  recall  the  three-dimensional  equation  for  the  conservation  of 
energy,  namely 

£ J #(*  /•  v*+«Vdv  = r (sV  + f*.  vVdv+f  ,(t.  v*-kV)da  , (7.16) 

dt  Jp*  ~ ~ Jp*  ~ ~ J^p4  ~ ~ 

* 

where  e is  the  internal  energy  density.  Suppose  in  addition  to  the  approxima- 
tion (6.6)  and  (6.7)  for  the  displacement  and  velocity  vectors,  we  adopt  the 
approximation 

K 


e"  = 8(e0,,t)+  £ (e3)N  ew(e“,t)  , e>o 

N=1  " 


(7-17) 


for  the  temperature  field.  Then,  for  a shell-like  region  bounded  by  the  surfaces 
(6.4)  and  (6.5),  from  the  energy  equation  (7*16)  we  can  derive  the  equation 
of  balance  of  energy  (3.12)  for  a Cosserat  surface  provided  we  make  the 
identification 

rP 


(*“  * * 1 3 

Xc  = pc  g"d0~ 


(7.18) 
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8.  Heat  flux  sectors  and  internal  energy. 

Suppose  the  ttiree- dimensional  shell-like  body  is  in  equilibrium  with 

* 

v =0  and  all  functions  are  independent  of  the  time.  As  a part  of  their  thermo- 
dynamic  restrictions  on  constitutive  equations,  Green  and  Naghdi  (1977)  have 
adopted  the  classical  inequality^ 


"3,  ' £ - 0 or  '£  ‘ £ - 0 » 


g = grad  0 


(8.1) 


for  all  time- independent  temperature  fields.  It  follows  from  (S.l)^  that 

(•Pi#  # q 

- j g^P  ■ g d0'  * 0 . ( 

or  ~ 

* 

With  the  approximation  (7-17)  for  0 and  the  help  of  (7.12)0  and  (7.13),  it 
is  seen  from  (8.2)  that 

K K 

E p^nVE-  S+  2 • §n  - 0 ( 

N=i  N N N=1 

for  all  equilibrium  displacement  and  temperature  fields.  With  the  above 
motivation,  we  add  the  inequality  (8.3)  for  all  equilibrium  states  to  the 
thermodynamic  inequality  (5*4)  which  was  derived  directly  from  two-dimensional 
postulates . 

Now  suppose  that  the  Cosserat  surface  C is  at  rest  with 


(8.2) 


(8.3) 


v = 0 , w = 0 


(8.4) 


for  all  time  and  with  the  deformation  gradient,  director  and  director  gradient 
each  constant  for  all  time.  Then,  by  (2.15) > P is  independent  of  t.  In  addition 
we  restrict  the  temperature  fields  to  be  spatially  homogeneous  so  that  9=  9(t), 

^Previously  (Green  and  Naghdi  1970)  in  the  case  of  an  elastic  shell,  an  inequality 
was  derived  using  (8.l)j_.  Because  of  the  different  thermodynamic  restrictions 
employed  here,  it  is  more  convenient  to  start  with  (8.l)0. 
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(8.6) 


9^  = 0 (N  = 1, . . . ,K) . Keeping  these  conditions  in  mind,  from  a combination  of 
i, 3 • 10 ) and  (3.1-3)  we  have 

K K 

p(r  + E r ) - div(q  + £ q ) = pc  . (8.5) 

N=1  N=l*'n 

In  view  of  (8.4),  no  mechanical  work  is  supplied  to  the  Cosserat  surface  C> 

Hence,  using  (8.5),  the  heat  supplied  to  a part  P of  the  surface  during  the 
time  interval  t^StSt0  is 

'I'o  f»  ^ ^ t 

M = I [ (r  + £ r )p  da  - (q  + E QjJ  • v ds]dt  = pc  c . (8.6) 

Jt1  JP  N=1  N JSP~  H-l**1  ~ JP  tx 

Suppose  the  shell  is  in  thermal  equilibrium  during  some  period  up  to  time  t 

with  constant  internal  energy  and  constant  temperature  Q.  We  assume  that 

whenever  heat  is  supplied  to  the  part  P of  the  shell  under  the  above  conditions, 

the  temperature  0(t)  throughout  the  part  will  be  increased,  i.e., 

t2 

[ 0].  > 0 whenever  W > 0 . (8.7) 

*1 

Provided  that  pc  is  continuous  and  remembering  that  p,  which  is  independent  of 
t,  is  positive  and  p is  arbitrary,  it  follows  from  (8.6)  and  (8.7)  that 


0(t)  - 0 > 0 whenever  «(t)  - > 0 


(8.8) 


for  all  t>t^. 


9.  Symmetries 

In  the  rest  of  this  paper  we  restrict  attention  to  the  situation  where 
K 1 in  the  earlier  part  of  the  paper,  and  for  convenience  we  set 


ei  ‘ ’ 


(9.1) 


The  inclusion  of  a second  temperature  field  0 allows  us  to  take  some  account  of 
temperature  variations  across  the  shell  thickness  but  the  more  general  situation 
in  which  K>1  can  be  dealt  with  in  a similar  way.  We  also  use  relative  kinematic 


measures  e ,y. ,K.  , instead  of  a ,d. , where 

CVp  1 Of  Ofp  1 lki 


eorP  ’ Vi 


d.  - D.  , K . = - „\. 

i l ta  to  iar 


(9.2) 


We  consider  the  form  of  the  Helmholtz  free  energy  function  in  (4.17)  which  is 
such  that  the  Cosserat  surface  models  the  main  features  of  a three-dimensional 
thin  elastic  shell  which  has  the  following  properties  in  its  reference  configura- 
tion: (i)  it  is  of  uniform  thickness  and  normals  to  the  middle  surface  meet 

the  major  surfaces  of  the  shell  at  points  which  are  equidistant  from  the  middle 
surface,  (ii)  it  possesses  isotropy  with  a center  of  symmetry  and  (iii)  it  is 
homogeneous  and  of  constant  temperature. 

We  choose  the  initial  director  D to  be  specified  by 


D = DA,  , 


where  D is  a nonzero  constant  in  line  with  (i)  above.  Then,  from  (4.11), 


A = - DB  , = 0 . 

arP  ap  ?a 

The  Helmholtz  function  (4.17)  may  now  be  replaced  by  the  different  function 


(9-3) 


(9.4) 


♦ = *(t'aP,VKicr’9’'>  J D,PIVV  ’ 


(9-5) 


which  by  virtue  of  (iii)  does  not  depend  explicitly  on  0°.  Moreover,  recalling 


35. 


(6.6),  our  Cosserat  surface  C models  a three-dimensional  body  which  in  its 
reference  state  has  a position  vector  R specified  by 


R*  = R + 9'^  DA.  . 

~ ~ -*0 


(0.6) 


At  this  stage  we  relax  the  condition  (ii)  slightly  and  assume  that  the  three- 
dimensional  body  is  transversely  isotropic  at  each  point  with  respect  to  the 
normals  A-  to  the  surface  9" = 0 (R  = R).  Then,  the  three-dimensional  energy 
function  is  form- invar iant  under  the  coordinate  transformation 

e3=-e3  , i°=9a(e1,92)  , 

where  9 any  other  set  of  three-dimensional  curvilinear  coordinates  with 


(9-7) 


9 on  the  middle  surface.  Corresponding  to  the  first  of  the  transformations  in 
(9-7),  we  assume  that  the  response  function  in  (9.5)  is  unaltered  in  form 


D-*-D  , d — - d , , 


the  other  variables  remaining  unchanged.  Thus 


(0.8) 


^ * ~^i 5 ia ’ ® 5 'D,‘DBa0,^y  ' 


= *(ec*’VSa’e’*;D>DVV  • 

Since  D/0,  ^ may  be  reduced  to  a different  functional  form  and  we  have 


(9.9) 


T <*0  i la’"’  ot$’~y 

if  redundant  elements  are  omitted.  To  complete  the  specification  that  the 

three-dimensional  shell  is  transversely  isotropic,  we  suppose  that  ^ is  an 

isotropic  function  with  a center  of  symmetry.  Moreover,  in  view  of  the  second 

^Complete  isotropy  is  considered  in  §10,  where  constitutive  coefficients  in 
the  linear  theory  are  identified. 


(Q.10) 


" 


part  of  the  condition  (i),  we  assume  that  the  major  surfaces  of  the  shell  in 
its  reference  configuration  arc  specified  by  9 - i } h in  (0.6',  where  h if  a 

constant  and  0'  = 0 is  its  middle  surface.  We  identify  the  .urfac«  of  c in  the 

« 

reference  configuration  by  K - R and  we  then  require  that  y is  ana tered  in  form 
under  the  transformation 


A3--A3  > --*3 


(9.11) 


The  transformation  i, 9 • 11)  implies  that 


D - - D , y0  --  , Y - Y » 

3 3 a a oor  00 


^aP  'aP  * ap  aP  ’ ^ aP  aP  ’ 


0 9 > • 


(9-12) 


Hence,  y in  (9.10)  is  subject  to  the  further  condition 


*(l'oP,DV3’'DVa’DK3a,'DKaP’9,"I>?  ; D ’'BaP’V 
= T(f'a3»DV3>DVa»DK3a,PKae,e,Dd  ; • 

After  imposing  the  restriction  (9. 13)  it  is  convenient  to  make  the  special 
choice  D=1  for  D which  implies  that  h is  the  thickness  of  the  shell.  Thus, 
(9.10)  reduces  to  the  different  form 


(9.13) 


* " ^eap,Y3’Ya,IC3a,K:ap,e^  : F«a’&/'  » 


aP 


(9.14) 


which  is  isotropic  with  a center  of  symmetry  and  which  is  also  subject  to  the 


restriction 


*(^P’Y3’"VSa’_'W0,''f  5 'W 
= *{VY3’VVW9’*5  VV  * 


(9.15) 


37. 


A similar  discussion  may  be  carried  out  for  the  entropy  flux  vectors  p ,p 

•%£)  '><1 


and  the  internal  production  of  entropies  £,»  , except  that  all  of  these  func- 
tions also  depend  on  d9/d9^,  d0/d9lJ'.  In  view  of  (9.3)  and  ( 9 • ) » the  entropy 
flux  vectors  may  be  reduced  to  functional  forms  of  the  type  exhibited  by  ^ in 
( 9 • ‘3 ) > including  dependence  also  on  dQ/dS11  and  d0/d9^*  Thus 


P®  = ,tc,  ,9,0  ; 


o ‘o'  a0”i’  ia 
o ~o 


de  d0 


; D,DB  ,A  ) , 


d^'  d0^  ' ^ 

. d9  d0 


(9-16) 


p = Pi^’V1^’0’*  ; — , ; d,db  ,a  ) . 

l 1 30  i ia  ^ d0M-  a0  ~Y 

Under  the  transformations  (9-8)  and  (9-12),  we  require  that 

p • g = p®^ea  , p^  • g,!  = p“a?/a9a  , §0  , ^0 


(9.17) 


be  unaltered  in  form.  After  setting  D=l,  it  follows  that  p®,p®  reduce  to  the 
different  forms 


p‘  = £ ; W ’ 

Pi  ~ Pi  (® fl>Yo>Y  ,K.  »K  ,0,0,  ■,  ^ > B ,A  ) , 

1 1 30  3 a 33  30  30 


(9-18) 


d0  d9K 

which  are  isotropic  with  a center  of  symmetry  and  which  are  also  subject  to  the 
restrictions 


VV'V  V'W9’'*’9*7881*’'8*799'*  : ‘ VV 
* P^-y3’VVV9’,,S9/a9*’*/a9“;  VV  ’ (9'1?> 

VV3’'V  V'‘<»e’9’'®’S9/S9VM/8^  ! 

■•5?(',«S'y3-yt,’,t3<.'K«9’e’<’99/99l‘’M/ie‘‘;  “ca'iy'  ’ (9-20> 

with  similar  forms  for  5,5^,  respectively.  Finally,  each  function  is  also  a 
hemihedral  isotropic  function  of  its  arguments. 


In  order  to  make  the  above  conditions  explicit,  we  limit  cur  attention  to 


1 


the  linear  theory  of  a thermoelastic  Cosserat  surface  which  is  unstressed  and 

t 


at  uniform  temperature  in  its  reference  configuration.  Then,  for  this  theory, 
the  position  vector  r and  the  director  d assume  the  forms 


r = R + u , d = DA_  ■»  6 , 


u = uM.  = u A*  , 6 = 5 = 6^A.  , 


(9.21) 


where  u,6  and  their  space  and  time  derivatives  are  small.  We  also  limit  our 

IV  iv  — — 

attention  to  a Cosserat  surface  which  models  the  main  features  of  a three- 
dimensional  thin  elastic  shell  ns  described  in  (i)-(iii)  above;  and,  hence,  we 
set  D=1  throughout  the  rest  of  the  discussion.  The  linearized  relative  kine- 
matic measures  are 

ecrP’VPia  ’ 


(9-22) 


where 


e _ = Mu  | . + uQi  ) - B „u-  = Mu  • A + u • A ) , 

c*P  a|P  p|cr  arp  3 ~>a  ~P  ~,  P ~a 


Pap  ' KaP + BapY3  ’ p3a  *3a " V\  * 


Vo  = 63  , Y„  = + u,  „ + b\iv  , = y^j 


(9.23) 


a 3> a a \ K3a  3j<* 


-p  = U„.  + b\  u + b\i  , + b\i  , -B  b\i  - Vi  , 

po8  3|oP  o|p  X a \|p  P aX  p 3 Ya|P 


3 


and  where  a vertical,  line  now  denotes  covariant  differentiation  with  respect  to 

the  reference  surface  using  Christoffel  symbols  calculated  from  A . 

crP 

We  assume  henceforth  that  the  reference  configuration  of  the  Cosserat 

surface  C coincides  with  its  initial  undeformed  configuration  and  thus  the 

vector  fields  representing  the  force,  the  director  force  and  the  internal 

director  force  are  all  zero  in  the  reference  configuration.  To  avoid  the 

introduction  of  additional  notations,  we  now  regard  the  force  vector  N and 

Mt  this  point  some  differences  will  occur  in  the  developments  if  we  wish  to 
include  the  effects  of  surface  tension  or  initial  stress  in  a membrane  theory 
of  a surface,  but  this  is  straightforward. 

39- 


the  director  force  vector  M as  infinitesimal  quantities  measured  per  unit  length 
of  curves  in  the  surface  S in  the  reference  conf iguration  of  c>  so  that  \>  is 
now  the  outward  unit  normal  to  curves  in  3 and 


v = \>  A'*  , N “>  , M = M** \j 

~ Qt~  ~ Qf  ~ ~ 01 


Moreover 


N0'  = iP'V  , M°  * , m = m1 

^ '•-'1  ~ r*s  1 ~ 


and  the  equations  of  motion  in  tensor  components  become 


N®*,  - pPv®+  pfP  = p(v  + kXw)  • AP  , 

|ar  o ~ ~ ~ 

vV+ V“Btpf3  * p<itl,1i)  -i?  - 


and 


M“K,  - V*+ pip  = p(k  v + kS*)  • AP  , 
I Of  ~ ~ ~ 


4°*,  -V^piP  '-1-* 

I Of 

M*3,  - V3  + pi3  = p(kXv  + k‘ w)  • A-  , 

| a ~ ~ ~i 


where 


(9.24) 


(9.25) 


(9-2o) 


(9.27) 


NP3  = yP  = mP  + bPm®3  , V3  = m3  - B M01®  , 

a arp 

f = f1A.  , l = A1 A. 
and  p denotes  the  reference  density. 

Turning  to  the  temperature  variables  we  replace  0 by  9+9,  where  ® ^ is 


(9.28) 


The  order  of  indices  N yrf*  has  been  used  by  Naghdi  (1972)  and  in  many  earlier 
papers  in  order  to  correspond  to  a usual  notation  in  shell  theory.  In  the  first 
paper  by  Green,  Naghdi  and  Wainwright  (1965)  the  opposite  order  was  used 

and_this  would  be  more  in  line  with  the  cooruinate  free  notations  of  the  form 
N N w = N10^  A,  sometimes  employed  in  the  current  literature. 

/v  «v  <y~\. 


constant,  and  assume  that  0 and  0 and  their  derivatives  are  small.  Then,  with 


P = P°A  , p = p®A  , 
~ ~cr  1 l~o 


(9.29) 


the  equations  of  balance  of  entropy  can  be  expressed  as 


pTl  = p(s  + 5)  - p“,  , 


(9.30) 


Also,  from  (3.3)  and  (3-9),  retaining  only  first  order  terms  we  have 


q=0p  , <1^  = 0 > r = 0s  , r1  = 0 . 


(9.31) 


We  recall  from  the  development  of  §4  that  most  of  the  constitutive  response 
functions  in  the  (nonlinear)  thermoelastic  theory  of  a Cosserat  surface  can  be 
expressed  in  terms  of  the  Helmholtz  free  energy  response  function  and  that  this 
is  mainly  due  to  the  use  of  the  energy  equation  as  an  identity  for  all  thermo- 
mechanical  processes.  A similar  situation  exists,  of  course,  in  the  linearized 
theory  and  leads  to 


* = *(ea0’VPia’e’*>  V£«) 


(9.32) 


K’0*  - i f(^-  * £-)  , ^ . v1 

SeB«  3pl« 

■n  - _ M - = 

ae  ’ \ d0  » 


P dY.  * 


(9.33) 


where 


N,0fp  = N'^*  = . 


(9.34) 


Further,  if  the  various  forces,  director  forces  and  the  entropies  all  vanish 
in  the  reference  configuration,  the  response  function  in  (9-26)  is  a 


quadratic  function  of  eap • This  function  is  to  be  a hemihedal 
isotropic  function  which  satisfies  the  condition  (9. 15)- 

ini 

In  order  to  obtain  an  explicit  form  for  \|i,  it  is  convenient  to  introduce 
the  invariant  surface  tensors 


E = e aAa8AP  , B = B aAa»AP  , 
~ ~ ~ ap~  ~ 

i ’ • £ - p3«p3» ■ 


(9-35) 


where  E,B,K  and  R are  symmetric.  Then, 


IV  /v 


* = *(e,j,k,p,r,y  ,0,0  ; B)  , 

* ini  ini  ivi  ini  ini  S ini 

where  ^ is  an  isotropic  scalar  invariant  function  of  its  arguments,  such  that 
it  is  a quadratic  form  in  the  variables  and 

V(E,-J,K,-P,R,v  ,9,-0  ; -B)  = t(E,J,K,P,R,v  ,6,0  5 B)  . 

INI  IV  IV  IV  |V  J IV  IV  IV  IV  IV  IN/  ^ V 

INI 

It  follows  (Rivlin  1955)  that  ^ is  a lineai'  function  of  the  following  39 
invariants 

Y?  , Vo  tr  E , (tr  E)2  , (tr  J)2  , tr  E2  , tr  J2  , 

J)  IN/  IN/  INI  INI  INI 

tr  JJT  , tr  K , tr  R , 


(9.36) 


(9-37) 


(9*38^ 


42. 


“ 


- 


0 , By,  , 0 tr  E , 

3 ~ 

2 

0 , 0 tr  J , 

Y-  tr  J B , v..  tr  E B tr  B , v,  tr  J tr  B , 

'3  ~ ~ 3 ~ ~ ~ 3 ~ ~ 

tr  E tr  J B , tr  E tr  J tr  B , tr  E tr  E B tr  B , 

***  ^ iv  /w  /v  »v  ^ m Arf 

(tr  J B)*"  , tr  J B tr  J tr  B , tr  E B tr  J B tr  B , 

******  >^*  ***  r**  ^**  ******  ******  /** 

tr  E B tr  J , (tr  E B)2  , tr  E J tr  B , 


* ***  i** 


tr  K B tr  B , tr  R B tr  B , tr  P tr  B , 

******  r**  ******  /**  ****** 


tr  EJB  , tr  BJJ  tr  B , 


***  <v  *** 


0 tr  J B , 9 tr  J tr  B , 0 tr  E B tr  B , 90  tr  B , 


0yo  tr  B , 0 tr  E tr  B , 0 tr  J B tr  B , 0 tr  E B , 

j ***  ******  ******  ^*  ****** 


(9.39) 


(9-*+o) 


(9.4D 


with  coefficients  which  depend  on 

(tr  B)2  , tr  B2  . 

***  *** 

Constitutive  equations  for  N,0(^,  (ani*  bence  may  be  obtained 

from  (9.33),  (9-34)  and  (9.38)  to  (9-41). 

It  remains  to  consider  the  constitutive  equations  for  PqjPj^ 
determining  these  we  must  also  use  the  identity  (4.8)  which  here  becomes 


(9.42) 


p(9+0)5+P0  = 0 > 


and  the  in.^quality  (8.3)  which  here  reduces  to 


(9.43) 


P*Sl  + £’&  + El*Si  *° 


43. 


(9-44) 


In  the  linearized  theory,  we  assume  that  p,^  and  5,^  are  linear  hemihedral 
isotropic  functions  of  degree  one  in  e .y^p^, 9,0,39/30^, 30/^  with 

coefficients  which  depend  on  A Q,B  such  that  p,p.  and  F,f.  are  subject  to 

orp  otp  **  ^ 

the  conditions  (9.19)  and  (9.20),  respectively.  If  we  retain  only  linear  terms 
in  (9-43^,  we  see  that^ 


§ = 0 


Moreover,  with  the  help  of  (9.44),  it  follows  that 


So  =-ao£-bo£Si  > gj.  =-ai£s-bi£ 


LSI  » 


= " ^2^  ’ P£\  = ' a3®br  E~  b3^  " c^tr  B E - d^  tr  J , 


where 


g = (ao/ae^A®  , g^  = (aa/ao^A®  , 


a0,.**,d3  are  scalar  functions  of  the  invariants  (9.42)  and 

b2  * 0 , aQ  * ° , bi  5 ° , aog  • g + (a^)  Bg  . & * 

for  all  g)§1. 


^Equation  (9-43)  then  determines  the  second  order  terms  in  5,  but  these  are 
not  required  here. 


(9-45) 


(9-46) 


(9-47) 


(9-48) 


V 
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10.  Thermoelastic  plate 

When  the  reference  surface  $ is  a plane  so  that  B=0,  the  results  of  the 
previous  section  simplify.  In  particular,  we  have 


* * “8a“^3«p3P  + 2"9A“Se«8v3 


* 2^V+2piA“Se«e9+2B2A“fi,W> 


+ B,62  * B432  , 


(10.1) 


and 


PTI  = - ( 0OV3  + PiA“0eap  + P3 0)  , p\  = - [ 02AafPpap  + P40 ] , 


<*P. 


E = -aoS  ’ El  =-bi§i  » 


(10.2) 


§ = 0 > P§x  = " b20  > P§x  = " b30  ' d3  tr  J > 


where  the  coefficients  in  (10. l) -(10.2)  are  constants. 

We  now  suppose  that  the  Cosserat  plate  characterized  by  (10. l)  and  (10.2) 

and  other  relevant  linear  constitutive  equations  models  the  small  deformation  of 

a (three-dimensional)  plate  of  initial  constant  thickness  h and  of  an  elastic 

material  with  the  following  properties:  Young's  modulus  E,  Poisson's  ratio  v, 

* 

conductivity  K , specific  heat  at  constant  volume  c,  uniform  density  p , coef- 

* _ * 
ficient  of  linear  expansion  a and  initial  constant  temperature  0.  Let  R be 

the  position  vector  of  any  point  of  the  (three-dimensional)  plate  in  its 

reference  state  given  by 


* 3 

R = R + eJA0  , 


~3 


(10.3) 


where  R is  the  position  vector  to  the  middle  plane  (0^  = O)of  the  plate,  and  let 
the  major  surfaces  of  the  plate  be  specified  by  0J = t 2 h • In  these  circumstances, 


45. 


I 


specific  values  have  been  suggested  for  most  of  the  constitutive  coefficients 
occurring  in  (10.1)  and  (10.2)  (see,  e.g.,  Naghdi  1972].  These  values  are 

“l  = a9  = 'l-2^  C » «2  = £'1-v)C  » “4  = (l-2v  C ’ 

n = viR  . rt  _ — ot  = /v  = ^ i l h 1 1 


(10.4) 


where 


<*5  = wB  , c*6  = = i(l-u)B  , <*3  = | tih  , erg  = ^ p.h 


r _ _Eh_  Eh  E 

" l-v2  ’ 'l2(l-v2)  ’ 


(10.5) 


^o  = Pi  = - » P2  = 


3 * 
Eha 

12(l-2v) 


(10.6) 


In  order  to  obtain  suitable  values  for  the  remaining  coefficients  we  make 
use  of  the  results  in  §7  where  all  quantities  now  refer  to  the  reference  state 


(10.3)  and  we  choose 


cr  = - 5 h , p = g h 


(10.7) 


For  a linear  elastic  solid  it  follows  from  previous  papers  of  Green  and  Naghdi 
(1970,1977)  that  if  only  linear  terms  are  retained,  then 

5*  = 0 , p*  = - (K/0)grad  9*  , 


* * * 
*_*  Ear  , . * . 0 c0 

p V = div  2,  + s~^~ 

9 


(10.8) 


where  u is  the  three-dimensional  displacement  vector.  We  use  the  formulae 
(7 -9) * (7.12),  (7.13) > (7-l4),  (10.8)  and  the  approximation  (7-17)  for  Q to 
obtain  the  following  values  for  the  coefficients  p^jP^d^ja^b-^b^b^: 


46. 
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V.'e  adopt  these  values  for  the  Cosserat  surface. 

To  complete  the  theory,  we  need  to  specify  values  for  the  assigned  force, 
the  assigned  director  force  and  the  assigned  supplies  of  entropy,  and  to  this 
end  we  consider  two  problems  below.  For  a discussion  of  values  of  the  assigned 
force  and  assigned  director  force  in  the  case  of  either  problem,  the  reader  is 
referred  to  Naghdi  (1972,  Ch.  E).  For  the  purpose  of  specifying  the  assigned 
entropy  supplies , we  suppose  in  the  first  problem  that  our  direct  theory  models 
a (three-dimensional)  plate  in  free  space  with  the  ambient  temperature  of  the 
surroundings  at  the  major  surfaces  of  the  plate  having  the  values 

e*  = at  03  - £h  + 0 , e"  = 0 at  0-'  =- Vh-0  . (10.10) 


The  temperature  across  the  faces  is  discontinuous  and  we  assume  a radiation 
type  surface  condition  of  the  form 


k*  = K(0  - 0+)  at  03  = £h  , 
k = K(0  - 0 ) at  0 = - yh  , 


(10.11) 


where  K is  a constant.  There  are  no  entropy  sources  in  the  plate  so  that  s 0. 

With  the  help  of  (6.8),  (7.10),  (7.11),  (10. 3),  (10.7)  and  (10.11),  we  find  that 

pso  = - K(20-  0+ - 0_)  , ps1  = - K[£hP0  - Vh(0+  - 0_)1  • (10.12) 


Collecting  together  values  of  ,pa,p®»  from  (10.2),  (10. o), 

(10.9)  and  110.12)  and  substituting  these  into  (9«30)  yields  the  following 
differential  equations  for  temperatures  0 and  i : 
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- 0 ^ (Yj  + A^e^)  + £-~]  +«hV20-H(20-  0+-  Q_)  = ° , (10. 13) 

0 

" * at  ['&{i-2v)  AOfepap  + + ^ • H[^  - ’Me+  - 9.)]  - Khd  = 0 , (10.14) 

10  9 


where  H=KQ  and  V"  is  the  two-dimensional  Iaplacian.  Apart  from  notation  these 
results  agree  with  those  obtained  from  a three-dimensional  approximation  by 
Green  and  Naghdi  (1970,1971)^  The  inequality  (9*48)  merely  yields  the  condition 
K £ 0 which  is  already  known  from  three-dimensional  theory.  If,  from  the  outset, 
we  had  restricted  the  discussion  to  a theory  with  one  temperature  field  0,  we 
would  only  have  the  one  temperature  equation  (10.13)  and  we  would  obtain  no 
information  about  the  temperature  variation  across  the  plate. 

In  the  second  problem  we  assume  that  the  region  - jh  5 9^5'h  is  a material 
interface  between  two  other  media  in  the  regions  9-  5 |h  and  9~  a^h.  The  tempera- 
ture and  heat  flux  (and  hence  also  the  entropy  flux)  are  continuous  across  the 
boundaries.  We  again  suppose  that  the  temperatures  in  the  surroundings  of  the 
interface  at  its  boundaries  have  the  values  (10.10)  and  we  now  suppose  that  the 


entropy  fluxes  at  these  boundaries  are 


k = k+  at  0"^  = ^h  + 0 , k = k at  9'  = - ^h  - 0 


(10.15) 


Because  of  the  continuity  of  temperature  and  in  view  of  the  representation 
(7-17)  in  the  interface  with  K = l,  we  choose 


9+  = 9 + £h0  , 9_  = 9 - ,Vh0 


(10.1b) 


Since  k is  continuous  at  9^  = ± £h,  it  follows  from  (6.8),  (7*9)>  (7*10),  (10. 3) 
and  (IO.15)  that 


Kxcept  for  a misprint  in  Elq.  (9-27)  in  197°  and  Bq.  (4.5)  in  1971.  The  value 
of  in  these  papers  should  read:  0-,  = - Eh3cr/l2  (l-2T|)  with  T1  being  Fbisson's 


ratio. 
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This  paper  is  concerned  with  thermomechanics  of  thir.  shells  by  a 
direct  approach  based  on  the  theory  of  a Cosserat  surface  comprising 
a two-dimensional  surface  and  a single  director  attached  to  every 
point  of  the  surface.  In  almost  all  previous  developments  of  the 
thermo-mechanical  theory  of  shells  by  direct  approach,  only  one 
temperature  field  has  been  admitted.  This  allows  for  the  character! 
zation  of  temperature  changes  along  some  reference  surface,  such  as  . 
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20  Abstract  (continued) 

the  middle  surface,  of  the  (three-dimensional)  shell-like  body,  but 
not  for  temperature  changes  along  the  shell  thickness.  A main 
purpose  of  the  present  study  is  to  incorporate  the  latter  effect 
into  the  theory;  and,  in  the  context  of  the  theory  of  s Cosserat 
surface,  this  is  achieved  by  a recent  approach  to  thermo:-. echani cf 
(.Green  and  Naghdi  1977)  which  provides  a natural  way  of  introducing 
two  (or  more)  temperature  fields  at  each  material  point  cf  'he 
surface.  t Apart  from  full  discuss  ion  of  theraomechanics  cf  shells 
and  thermodynamicalVrestrictions  arising  from  the  second  law  of 
thermodynamics  for  shells , •''attention  is  given  to  a discussion  of 
symmetries  (including  material  symmetries)  ana  thermal  effects  in 
the  nonlinear  theory  of  elastic  shells  with  detailed  discussion 
of  the  linear  theory  of  elastic  plates.  — 


